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Abstract 

New techniques for evaluating the closed time path efFecive action for non- 
equilibrium quantum fields are presented. A derivative expansion is performed 
using a proper time kernel. Applications relevant to the scalar field theory 
of thermal inflation are discussed and dissipation terms resummed. The ef- 
fective action of the electromagnetic field is also considered. In this case the 
leading term can be related to the conductivity of a plasma and has a simple 
interpretation in terms of the classical Drude theory of conductivity. 

Pacs numbers: 05.70.Ln, 98.80. Cq, ll.lO.Wx 
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Typeset using REVT^ 



I. INTRODUCTION 



Many situations of interest in Particle Physics and Cosmology are close to thermal equi- 
librium with fields which have slowly varying ensemble averages. Two examples in the early 
universe are where the time evolution is due to the expansion of the universe or the slow 
change of an infiaton field. In is important in such situations to understand both how dis- 
sipation affects the dynamics and how the how the dynamics, in turn, reacts back on the 
thermal ensemble. 

In quantum field theory, the behaviour of ensemble averages can be usefully encoded in 
an effective action. For systems in thermal equilibrium an imaginary time formulation is 
usually used but this is only appropriate for slowly varying spatial fields @,|]. For 

systems that are not in thermal equilibrium, the effective action can be introduced within 
the closed time path formalism invented by Schwinger which allows ensemble averages to 
evolve in time The relativistic formulation |]^,^ allows a full account of the dynamics 

of nonlinear quantum fields and thermal fiuctuations. 

The case of slowly varying ensemble averages lends itself to a derivative expansion of 
the effective action. Terms that are linear in time derivatives are responsible for dissipative 
effects, where the short wavelength modes behave as a thermal bath damping the long 
wavelength average. Early work was done on this problem by Hosoya et al. 0, who were 
studying the problem of reheating in the inflationary universe. The flrst derivative terms in 
the closed time path formalism have been studied for a self-interacting scalar fleld by Gleiser 
et al [IlOi . Higher order derivative terms have been obtained by Das and Hot , although 



they restricted themselves to terms quadratic in the flelds. Berera et al. [|12l have further 
examined the flrst derivative terms and done some preliminary work on higher derivative 
terms. Their work has been largely stimulated by the warm inflationary scenario, where the 
evolution of the inflaton fleld is damped by thermal interactions |T^-|T5|. 

The flrst problem in a derivative expansion is that free fleld theory is not a good place 
to begin. This is because a non-interacting system has no way to re-establish thermal 
equilibrium after an external disturbance, and so successive disturbances will drive the 
system further and further from equilibrium. This problem could be avoided by considering 
only a small time interval [|T6|, but the prefered method is to begin with an equilibrium theory 
with a dressed propagator which damps down the effects of the external forces |T^,|T^JT^ . A 
new justiflcation of this approach based on the effective action will be given in appendix A. 

We shall see that the propagator has a simple representation in terms of a proper time 
kernel. Unlike the zero temperature case, the kernel is nonvanishing in both forwards and 
backwards proper time directions. The use of the proper time kernel simplifles the derivative 
expansion of the effective action. 

Application of the derivative expansion to a scalar fleld coupled to a heat bath conflrms 
that the derivative expansion is an expansion in powers of 7~^5t, where 7 characterises the 
damping effects of the particles in the heat bath. Resummation of the derivative expansion is 
possible in the case that the second and higher derivatives of the averaged fleld vanish. The 
results indicate that the resummation introduces logarithmic terms in the high temperature 
limit . 

In an attempt to understand the necessity for including damping terms in the derivative 
expansion the formalism is applied to Quantum Electrodynamics in section V. Dissipative ef- 
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fects at non-zero temperatures are well known physical phenomena. Constant electric fields, 
for example, drive currents through the plasma. These can be calculated from the nonequi- 
librium effective action. The current can also be estimated by elementary Drude theory and 
demonstrates how the damping is related to collisions between accellerated electrons and 
thermal electrons or photons. 

The quantum field theory conventions used are those of Weinberg ||19[ , in which h = c = 
1, the metric signature is — h ++ and the Feynman propagator is (fc^ + — ie)^^. 



II. CLOSED TIME PATHS 

The closed time path formalism provides a definition of Green functions and the effective 
action of a nonequilibrium quantum system. The path integral version used here follows 
Calzetta and Hu p. 

At the initial time ti the probability distribution of the quantum field is given by a 
density matrix p[0i,02]- The generating functional is then defined by a path integral over 
the two fields 0i and 02 which are given the same values at the final time t/, 

Z[Ji,J2] = I rf/i[0i] / rf/i[02]p[0i,02]e*^[<^^]+^/-'^*^-^^*[^^l-^/-'^*^ (1) 



where 5'[0] is the classical action and two sources Ji and J2 have been introduced. The 
generating functional has an important symmetry, 

Z[Ji,J2] = Z[J2,Jl]* (2) 

and, with suitable normalisation, Z[J, J] = 1. 

The generating functional is practically the same as any generating functional for two 
scalar fields except for the special boundary conditions. This similarity can be made more 
explicit by introducing a metric Cab = diag(l, —1), then .P = — J2. The effective action is 
therefore 

m = -tiog z-jy^ (3) 

where contracted indices a indicate an integral over x |20| and J"^ is the solution to 

(4) 



As usual, 

5T 



-r (5) 

The ensemble average (pc is defined by 0c = 0i = 02, which occurs when Ji = J2. 

The effective action can be approximated in various ways. A modified version of the one 
loop approximation will be used here, where T ^ S + W and 

= ilogdet(l + SGD). (6) 
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S contains dependent mass insertions and Gd is the equilibrium (J^ = 0) propagator. 
Specifically, 



(7) 



S(t)a6(t)b 

where Sj contains the interaction terms in the action and (in momentum space) 

{Go')\={e + ml)6\ + Eo\ (8) 

where Sq is the sum of equilibrium self-energy corrections to the propagator. 

Using dressed propagators avoids divergent terms arising in the derivative expansion 
TO| , P^ . In particular, it is the dissipation introduced by the mass correction Sq which bal- 



ances the effect of the external sources. A derivation of the corrected one loop approximation 
(I) is given in appendix A. 

Some general statements can be made about the form of the equilibrium propagator G. 
It will be useful to review these properties of G in order to introduce some matrix notation 
which will be used in subsequent sections. 

The propagator can be obtained from the generating functional by 



6J^{x)6J\x') 



(9) 



Ja = 



The symmetry (0) of the generating functional implies (for bosons) 



Gnix, x') = Gnix', x) = -^22(0;, x')* = -^22(0;', x)* (10) 

Gi2{x,x') = G2i{x',x) = -Gi2{x',x)*. (11) 

Furthermore, there are two physical conditions. Causality implies (see [^), 

Gil + ^22 + ^12 + G21 = (12) 

Thermal equilibrium imposes periodicity in the imaginary time direction, 

6-21 (x, t, x', t') = Gi2(x, t + i(3, x', t') (13) 

Fermion propagators are antiperiodic. 

The symmetry and the causality conditions reduce the number of independent functions 
needed to describe the propagator. For example, it is possible to write the propagator in 
terms of functions G^{x,x') and G^{x,x') [|T^], 

Gn{x,x') = +e{t-t')G>{x,x') +e{t' -t)G<{x,x') (14) 

G'i2(x, x') = -e{t - t') G>(x, x'Y + e{t' - t) G<(x, x') (15) 

G2i{x,x') = +e{t-t')G>{x,x') -e(t' -t)G<{x,x')* (16) 

^22(0;, x') = -e{t - 1') G>{x, x'Y - e{t' - 1) G'<(x, x'Y (i7) 

where G^{x^x') = G^{x',x). 
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In thermal equilibrium, the propagator can also be expressed in terms of the spectral 
function 

p = tG2i - tGi2- (18) 
In momentum space, the symmetries of the propagator imply that it takes the form 

G\ = la6\ + ^N%ep (19) 
where a = Gu — G22 and e(co') is the sign function. The periodicity in imaginary time gives 

with the upper signs for bosons and the lower signs for fermions. Two important properties 
of the matrix are A^^ = 1 ^^^1 ]V"^(_lj) = Nb^'^u). 

The general form of the propagator can now be compared with the dressed propagator 
(P). We see immediately that the self energy must have the matrix form 

T.^\ = 5m\k)5\-i^{k)N\ (21) 

where 5m? and 7 are real functions. If we set m? = + 5m?, then 

= (A;' + n?)5\ - i^N\ (22) 

Imposing G^Gd = 1 enables us to deduce that 

p = i^k"^ + + ie7)~^ - i(A;^ + - ^€7)"^ (23) 
a = {k'^ + n? - i-fY^ + {k'^ + m^ + i-iY^ (24) 

The poles of the spectral function p occur when uj = ±ujk ± iT, where ujk = (|kp + m^)^/^ 
and 

7(A;) = 2uje{uj)T{k) (25) 

The function r(A;) can be interpreted as a decay width. 

The propagator is often written in a hybrid form which depends on spatial momentum 
and time |[10|. The inverse transform can be evaluated by using equations ([T9|), ([20| ) and 
(H), then for t > 0, 

G>(k, t) = -^ (fiuk + 2r)e*(-'=+^r)* + (1 + fiuk - ^r))e-^(-'=-^r)*) (26) 

This is a very unwieldy expression to integrate, especially when there two or more propaga- 
tors involved in a perturbative diagram. An alternative method for evaluating integrals of 
propagators is given in the next section. 
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III. THE PROPER TIME KERNEL 



A new idea will now be introduced to simplify manipulations involving the propagator. 
This is the introduction of a proper time kernel, related to the kernel which was introduced 



into zero temperature field theory by Schwinger It is also related to the use Feymann 
parameters for the evaluation of Feynman diagrams. 

The proper time kernel Kab{x, x', r) is a 2 x 2 matrix, although the matrix indices will not 
be written out explicitly. For convenience in switching between configuration and momentum 
space the bra and ket notation will be used and the kernel regarded as an operator K{t). 
It satisfies a Schrodinger equation, 

idrK = AK + i6{T) (27) 

where AG = 1. The proper time kernel generates the propagator by 

G{x,x') =% \ dr 1^x'\K{t)\x) (28) 



oo 



Unlike at zero temperature, the proper time integral covers past and future values. 

For the dressed propagator G^i, the momentum space version of A = is given by 



equation (|2^) . There is a unique kernel which satisfies (^) and for which the integral (pSD 
converges. 



^D = H^ + <^))e ''^ (29) 



i^D(fc,r) = Ae-'(^'+'"')" (30) 



A = i(iV + e(r))e-^l^l (31) 



In momentum space. 



where, 



The matrix was given by equation (pO|) in terms of the distribution function /. 

It is instructive to examine the behaviour of the kernel at the inverse temperature ^ ^ oo. 
One finds 

^""(^^'"o'' 0(-r)e^O' ^^^^ 



where Q{t) is a step fuction. When substituted into ( PBD one recovers the Feynman and 
Dyson propagators at zero temperature, with a ^7 in the places where you would normally 
place the ie. 

The non-eqilibrium proper time kernel for the operator A = A/j + S can be obtained 
by perturbation theory. The usual expressions involving time ordered products cannot be 
used because of the acausal nature of the propagator. Instead, it is helpful to define the 
convolution product for functions of r, 

dr'f{T-r')g{r'). (33) 

-00 
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The kernel satisfies an integral equation 

K = Ko-iKoo EK. (34) 

The solution to the integral equation is 

K = KDo{l + iJ:KD)-^ (35) 

where repeated use of the convolution product is used in the expansion of the inverse. The 
nonequilibrium propagator can be recovered from (^), which allows us to write the corrected 
one loop effective action (^ as 

W = - rfr tr log (1 + iEKD) • (36) 

Z J— oo 

If the logarithm is expanded it is possible to use the composition property Kd o Kd = tKd 
and the cyclic symmetries of convolutions to write this in a form resembling Schwinger's 
zero temperature result, 

i 

W = - dTT'hiK + constant. (37) 
The one loop divergences show up at small values of r, but the proper time approach is 



especially well suited to trouble-free analytic regularisation schemes p2 



To derive the derivative expansion of the effective action we suppose that = 
S(x)5(x — x'). The full expression for a term (EKo)^ in a momentum space basis is 

{k'HEKorik) = {-iri:{idu)KDO...oi:{idk)KD5{k-k') (38) 

where Kr, = K£,{k,T) on the right hand side. In configuration space, 

(x|(Sir^)"|x) = (-z)" j + idk)KD o . . . o S(x + tdkWo (39) 

The derivative expansion is obtained from a Taylor expasion about S(a;). However, we first 
introduce the expression (^) for K^. If we take w? to be independent of k it is possible to 
move the exponential factors around to leave, 

{x\{J:KdT\x) = {-lY j ^0je-''"'"-''='"/2(S(x + 5)v4o...oS(x + 5)A)e-''='"/' (40) 

where the 5 in the j'th position along the convolution product is the operator 

5{Tj) = idk + 2kTj - kr (41) 



It is convenient to adopt the notation p6 



d'^k 



(...)= / -U-tre"'^'"/' . . . e-'^'"/^ (42) 



(27r)4 

then the one loop action can be written in a compact form 
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W = ^ fd^x r dre-'""'^ {log (1 + + 6) A)) (43) 

Z J J — oo 

where the terms expand to the full expressions (0). 

A useful property of the one loop action is that the mass can be shifted by an arbirary 
function X(x), 



W 



- f d*x r dre-'('"'+^)"(log (1 + z(S(x + 5) - X)A)). (44) 

2 J J —oo 



This follows from the properties of the proper time kernel. 

The effective field equations are obtained by differentiating the effective action. The one 
loop contribution is therefore 

SW 1 , __2_9S 



- r dTe-'^'^{^A (1 + + 6)A)-^) (45) 

/ J-oo 0(pa 

With the physical condition = 0i = 02, the function S is related to the interaction 
Lagrangian Ci{(f)), 

and we use the shifting property (^) to get the main result of this section 

where = m? + X and = diag(l, —1). 

The derivative expansion of the effective field equations is obtained by expanding in 
powers of 5 about x and evaluating the resulting matrix elements. The matrix A was given 
in equation (|3TD. 



IV. SCALAR FIELD THEORY 

Now we turn to an explicit calculation based on a semiclassical scalar field interacting 
with a heat bath of weakly interacting scalar particles. The semiclassical field will be re- 
garded as varying slowly in time so as not to upset the thermal equilibrium of the heat bath. 
This is an idealisation of the situation encountered in the critically damped evolution of an 
inflaton field. 

The choice of interaction Lagrangian for the semiclassical field and the heat bath fields 
7] will be 

C, = Ig'mv' (48) 

The self-interaction terms of both and rj will be assumed to be of order g'^, which means 
that the one loop effective action comes predominantly from the rj particles. 
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The dressed propagator in models of this type has been studied by Berera et al. |T2[ . A 
reasonable approximation can be made by evaluating the self energy at zero momentum. In 
the high temperature limit T ^ m, for example, 

m^^ml + (49) 

7 « ^ (50) 
' 647r ^ ^ 

for the theory with $ = 0^. 

The one loop contribution to the effective field equations reads 

The shifting property of the effective action allows us to replace by M'^{(j)) = m^ + (yf^$(0). 

We shall try to recover all of the terms invoving $ = 9^$ and neglect and higher 
derivatives. Thus, since $,^^5'^ = — 

where, 

oo „ 

dre-'^' ^z"(la3A(5^A)") (53) 

-oo 

The evaluation of the matrix element and proper time integral are described in appendix B. 
To leading order in 7, 

1 d:f ^^^^ 



(27r)3 2uJk (27)- 

where ujI = |kp + M^. The term with coeficient ci agrees with the results of Berera et al 



g. 

The behaviour of the sum can be understood a little better when it is expressed as an 
integral. To begin with, let $ be negative. Set 

The one loop term becomes 

(^)^=/$7_°^n(/.)e--/(^^*)ci;. (56) 

Analytic continuation can be used for large positive $, with a pole at $ = 2'yT/g'^. 

The integral can be further simplified in the high temperature limit T ^ M. Using 
standard techniques for high temperature expansions gives 
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for negative /i. Inserting this into ( ^61) gives 
(5W\ g^<l>'T^ g^^ 



M/^)~-^log(l-e'^'^) (57) 




where is the digamma function. For small $, the result gives 



V(l) (58) 




167r2 27T ^ \2-fT 



log . (59) 




For comparison, the first term in the series expansion (|52D gives a similar result but contains 
no logarithm of $. Beyond the one loop approximation, the two-particle irreducible diagrams 
are of order g^, but the presence of a leading logarithm indicates that resummation of the 
one loop approximation gives the largest contribution to the effective action can be used 
consistently. 

Higher derivative terms can be obtained by similar means. The simplest term with two 
derivatives would be 

= --g'<^' / rfre-^"^ ^|(ia3AWA)X^, (60) 
Z J —00 

Using the methods described in appendix B, this reduces to 

U^A 24 

where Zi = —Ii[2ujdf] and Z2 = Isl^f + 1] in terms of the integrals defined in appendix B. 
The second term agrees with earlier results for spatial derivatives in reference 0, whilst for 
small 7, 

with / = f{uJk)- The first term is the dominant time-derivative term for small 7. There are 
additional contributions from two-particle irreducible diagrams at order g"^ according to the 
number of vertices, but these diagrams give fewer inverse powers of 7 and do not affect the 
leading term. 



Berera et al. |T2[ have claimed that the leading term is absent if the homogeneous limit 



(k — > 0) is taken after letting F — 0, as suggested by Gross at al. |2j. However, they also 
find that the first derivative terms diverge when the calculation is done in that order. In 
the method used here, all of the terms are recovered even when the fields are not homoge- 
neous. In the next section, the situation in quantum electromagnatism is investigated to 
help understand first derivative terms. 
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V. QUANTUM ELECTROMAGNETISM 



The effective action for the electromagnetic field at zero temperature contains corrections 
from vacuum polarization. The proper time kernel for a constant field can be obtained 
exactly the first correction to the effective action being the Euler Heisenberg term. 

Dissipative effects prevent such a simple treatment of the nonzero temperature problem. 
The leading term is due to the current induced by the electric field, which will be calculated 
here from the derivative expansion. The result will be compared to the electrical conductivity 
calculated by elementary transport theory p8 . 



The proper time kernel for electrons in a background elecromagnetic field F^'^ satisfies a 
Schrodinger equation with the operator 

A = -{d + teAy + m^-^ea^^F^"' (63) 

where 0"^,^ = |«[7^,7jy]. Comparison with the free theory shows that 

S = 2eA^k^ + e^A^A^ - \ea^,F^^ (64) 

In thermal equilibrium, the electron propagator contains mass corrections of order e^, 
but the first disipative terms arise at order e^. These have been discussed in detail by Le 



Bellac 29 



The corrected one loop contribution to the effective field equation is obtained from (^5|), 
/ dTe-''^^{\a^{2ek^ + iea^,k'')A{l + iT.{x + 6)A)-^) (65) 



5A^' 2 J- 
In this case, since 

5V 



-J,, (66) 



SAi' 

the one loop term can be viewed as an induced electric current. 

The calculations can be simplified by using Schwinger p5| gauge in which 

^^ll,^l2■■■^^n — ^^"qJY-^/i2Ml,/i3---Mn (67) 

To leading order, 

E = eF^.d^k'' - la^.F'^" (68) 
The leading term can be evaluated using the formula given in appendix B. The result is 

5W , , 

where the electric field is E = —F^^ and, for small constant damping 7, the conductivity 

r d^k 1 



a 



7 J (27r)3 e'^'^fc + 1 



(70) 
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with Lul = |kp + 771^. 

The result for a is of order e~^. This can be understood in terms of elementary Drude 
theory. In a plasma of density n, if an electron accellerates for a time r between collisions due 
to the applied electric field and then moves off in a random direction, the Drude conductivity 
will be 



ne^r/m. (71) 



The relaxation time r can be infered from the location of the poles of the dressed propagator 
uJk + ^r. According to equation (^), r = F^^ = 2m/7, which gives a result in agreement 
with (|70|). In terms of the collision cross section ar, the relaxation time is roughly l/{nv(7T)- 



Since the cross section is of order e^, the value of a is of order e ^. 



VI. CONCLUSION 

We have seen how the derivative expansion of the nonequilibrium effective action can be 
used to approximate the dynamics of a quantum field interacting with thermal radiation. 
The importance of damping in the thermal medium is very apparent and we have seen that 
the reason for this can be understood in terms of scattering theory. 

When the damping is small, the dominant time derivative terms at order n are always 
associated with n inverse factors of the damping coefficient 7. This implies that successive 
terms in the derivative expansion will be decreasing in size if the n'th derivatives are of 
order 7" or less. Exceptions can occur at phase transitions where infra-red problems affect 
the momentum integrals. Furthermore, resummations of the expansion series can lead to 
logarithmic terms, as we saw for the power series in 7"^$ in the high temperature limit. 

The derivative expansion for external electromagnetic fields in Quantum Electrodynamics 
which was discussed in section V can easily be extended to quadratic order where the 
results can be compared to zero temperature results obtained from the Euler-Heisenberg 
Lagrangian. The derivative expansion will give non-linear corrections to the conductivity in 
this case. It should also include the Hall effect. 

The proper time kernel was useful in helping simplify the derivative expansions but it 
can also be applied in other contexts. It provides a natural way to extend the Feynman 
parameterisation of the propagator to non-equilibrium Feynman diagrams which could then 
be used to evaluate multi-loop diagrams. It would also be valuable to have the exact proper 
time kernel in the presence of a constant background electromagnetic field, since it is known 
at zero temperature Unfortunatly, this appears to be a highly non-trivial problem. 
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APPENDIX A: THE ONE LOOP EFFECTIVE ACTION 



The one loop effective action can be derived from tlie generalised effective action of 
Calzetta and Hu [§ (also developed in refs. p3|j2^ ). This form of the action uses a source 
Jab quadratically coupled to the field (pa to generate the propagator Gab- The effective action 
r[0a,Gbc] has the perturbative expansion 

r[0„ Gbc] = S[ct>a] - \ logdet(G,fc) + l^^'^'Gab + r2 (Al) 

where r2 is the sum of two-particle irreducible diagrams (i.e. they cannot be disconnected 
by removing two lines) in a theory with 'shifted' vertices and 

0(paO(Pb 

The effective action r[0a] can be obtained by replacing the propagator with Gbc = Cfeci^a], 
obtained by solving 

(A3) 



?>G 



ab 



which implies 



{G-'y = A^' + (A4) 

(>Gab 



Hence 



r[0j = S[(Pa] - \ logdet(G,,) + E (A5) 

where Gab = Gab[4'a] and E is a two-particle irreducible remainder term. 

The dressed equilibrium propagator Gd is defined to be the value of Gab when cj) = 0, 
hence 

where S = A[0]— A[0] and E' is another remainder term consisting of one-particle irreducible 
diagrams which vanish when = 0. The corrected one loop approximation consists of 
dropping E and E' , 

r[0j ^ S[<P\ + logdet(l + IIGd) (AT) 

The error terms have at least two shifted vertices. In theory, they would be of order A^0^ 
at most, whereas the corrected one loop term is still valid for Xcj)^ large. Some examples of 
graphs contributing and not contributing to the one loop approximation are shown in figure 
1. 
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FIGURES 





(a) 



(b) 



FIG. 1. Feynman diagrams for a scalar field with interactions, (a) shows an example 
of a diagram contributing to the dressed one loop approximation used in the text; (b) shows a 
two-particle irreducible diagram with shifted vertices not included at one loop. 

APPENDIX B: USEFUL FORMULAE 

The derivative expansion produces matrix elements containing products of convolutions, 
derivatives and matrices. These can be evaluated with the methods described below. To 
begin with, consider the proper time convolutions. The convolution product is both asso- 
ciative and cyclic. We can usefully define operators d and d'^ acting on functions of proper 
time 



df{n) =Tif{Ti), f o d^g = (df) o g. 



fBll 



The operator d satisfies d{f o g) = {df) o g + f o dg. 

Derivatives with respect to momentum appear in the operators 6 (0). These operators 
can be expressed in terms of d and d^ by 



S=- (d^D - DU 



(B2) 



where 



D = dk + ikr, 



-dh + ikr. 



(B3) 



Note also that k = —^i{D + D'^)/t. The operators D and behave like anihilation and 
creation operators inside the matrix elements (|^), and satisfy 



(B4) 



The strategy for evaluating the matrix elements involves moving the D operators to the 
right and the operators to the left. 
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The proper time kernel depended on a matrix A defined in equation (pTI). Convolutions 
of A can all be obtained from the formula 

T) ^ Til ' 

d'^A o rA = (B5) 

(n + m + 1)! ^ ^ 

As a simplification we shall take the damping function 7 to be constant. This can be justified 
as long as we are interested in only the leading term in the limit 7^0. (Derivatives of 
7 can be retained if required.) Derivatives of A can be then be dealt with by using the 
formula, 

A- o (d^^A) o A" = f ^^^-^^L SnUNd)A - SUNd)A j (B6) 



where 



o r \ ^ f n\ in + m - r)\ ,„ 

Snm{x) = E J ^ r-r^ 27a: B7 

^ \ r J n\m\ 



r=0 



These formulae reduce the matrix products to expressions of the form d^Nd"'A and 
d"'A dl^N. The symmetry under uj —uj implies 

{a^{dlN)N'P+^S{Nd)A) = -{asS{Nd)A{d;;^N)NP+^) (B8) 

where S{x) is any polynomial. After substituting for the matrices A and A^, 

{as{dPN)NP-''SiNd)A) = 2(-ir((9£/)eP+i5(-9,)^^+) (B9) 

where 6'+(r) = e"'^!'^'. 

Now consider the proper time integrals contained, for example, in equation (|47|). These 
produce integrals of the form 



drizry-'e-'"^ ^ aas{d:N)N-^' S{Nd)A) = (BIO) 



where 



m = dri^rr-' J ^e-(^^+™^)^^4^)M^)- (Bll) 
The integral Is can be evaluated by noticing that 

dre-*('''+'"')"^+(r) = e{uj)p{k) (B12) 



where p is the spectral function (pSf ). For s > 2, Ig can be obtained by differentiation with 
respect to m^. For s < 2, analytic continuation must be used to obtain, 

^.l''l-fJ^^/|^l''P'-M-)pWM-) (B13) 
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This result is valid for arbitrary 7. A useful property of the result is that Is[duih] 
In the limit 7^0, the spectral function reduces to a delta function and 



r(|-s)i (27r)3' ' 2uJk 



with ujI = |kp + m^. 

As an example, consider the evaluation of 

M„ = (|(T3Ao(5,A)") (B15) 

The S can be expressed in terms of D using equation ([B2D. When moving the D's to the 
right, they can contract with other D's or act on the A matrices. For small 7, the leading 
contribution comes when all of the D's act on a single A matrix. After using equation ( |B5p , 

M„ = ±{-ir-^ (f )" ^^^'y^^-.f (iaaA^^- o (d^A) o (B16) 

Now equation (p6|) can be used. The following summation is needed, 

E(-l)^ ^ 1 ' ^2(n-,)2.(x) = (47x)"^„(x) (B17) 



j=0 



where 



Hence, 



n 1 

Sn{x) = j:-M' (B18) 

T- = 



M„ = i2t^r''aasidZN)N''+'Sn{Nd)A) (B19) 



The proper time integral of M„ used to define the coefficients c„ (|53D can be read off from 
equation ( |Bl(j| ), which reduces to c„ = (— in the limit 7 — > 0. 



16 



REFERENCES 



L Dolan and R Jackiw, Phys. Rev. D 9, 3320 (1974) 
S Weinberg, Phys. Rev. D 9, 3357 (1974) 

I G Moss, D J Toms and A Wright, Phys. Rev. D 46 1671 (1992) 

I G Moss and S J Poletti, Phys. Rev. D 47, 5477 (1993) 

J Schwinger, J Math Phys 2, 407 (1961) 

L V Keldysh, Zh. Eksp. Teor. Fiz. 47, 1515 (1964) 

G Z Zhou, Z B Su, B L Hao and L Yu, Phys. Rep. 118, 1 (1985) 

E Calzeta and B L Hu, Phys. Rev D 37, 2878 (1988) 

A Hosoya and M Sakagami, Phys. Rev D 29, 2228 (1984) 

M Gleiser and R O Ramos, Phys Rev. D 50, 2441 (1994) 

A Das and M Hott, Phys. Rev. D 50 6655 (1994) 

A Berara, M Gleiser and R O Ramos, Phys Rev D 58, 123508 (1998) 

I G Moss, Phys. Letts. B 154, 120 (1985) 

A Berera, Phys. Rev. Lett. 75, 3218 (1995) 

A Berera, Phys Rev D 54, 2519 (1996) 

M Asprouh and V G Gonzalez, Proceedings of the 5th International Workshop on Ther- 
mal Field Theory RegenshvLig, Germany, (1998), |hep-ph/9811469 . 



E Wang, U Heinz and X Zang, Phys. Rev. D53, 5978 (1996) 
P R Parwani, Phys. Rev. D45, 4695 (1992) 

S Weinberg, The Quantum Theory of Fields, Vol. 1 (Cambridge University Press, Cam- 
bridge UK 1995) 

B. S. DeWitt, Dynamical Theory of Groups and Fields (New York: Gordon and Breach 
1965) 

J S Schwinger, Phys. Rev. D 82 664 (1951) 

J S Dowker and R Critchley, Phys. Rev D 13 3224 (1976) 

I D Lawrie, Phys. Rev. D 40, 3330 (1989) 

G Amelino-Cameha and S-Y Pi, Phys. Rev. D 47 2356 (1993) 

J S Schwinger, Sources, Particles and Fields (Addison- Wesley, New York 1973) 

I. G. Moss and W. Naylor, Class. Quantum Grav. 16 (1999) 2611 

D J Gross, R D Pisarski and L G Yaffe, Rev. Mod. Phys. 53 43 (1981) 

C W Misner and D H Sharp, Phys Letts 15 279 (1965) 

M Le Bellac, Thermal Field Theory (Cambridge University Press, Cambridge UK 1996) 



17 



